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We present four estimators of the shared information (or interdepency) in ground-states given 
that the coefficients appearing in the wavefunction are all real nonnegative numbers and therefore 
can be interpreted as probabilities of configurations. Such ground-states of hermitian and non- 
hermitian Hamiltonians can be given, for example, by superpositions of valence bond states which 
can describe equilibrium but also stationary states of stochastic models. We consider in detail the 
last case, the system being a classical not a quantum one. Using analytical and numerical methods 
we compare the values of the estimators in the directed polymer and the raise and peel models which 
have massive, conformal invariant and non-conformal invariant massless phases. We show that like 
in the case of the quantum problem, the estimators verify the area law with logarithmic corrections 
when phase transitions take place. 
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As is well known, in quantum mechanics, for a pure 
state, if C = A + B is a bipartition, the von Neumann 
entanglement entropy S q is defined as 



S q (A) = S q (B) = -Tv(p A lnp A ) 



(1) 



where p A = Trg(p) and p is the density matrix related 
to the ground-state wave function. For one-dimensional 
spin systems defined by Hermitian Hamiltonians, if the 
lengths L and I of C respectively A, are large one has the 
area law S q stays finite if the correlation length is 
finite. If the system is gapless and conformal invariant, 
one gets logarithmic corrections 

S q {l,L)^ 1 \nl + C (L»l), (2) 

and the finite-size scaling behavior 

S q {l,L) = 7ln(Lsin(7r//L)/7r) + C, (3) 

where 7 = c/6 for an open system and 7 = 2c/6 for 
periodic boundary conditions (c is the central charge of 
the Virasoro algebra) Q. The factor 2 in the latter case 
appears because the systems A and B have two common 
boundaries. C is a non-universal constant. The relations 
©-© have been checked analytically and numerically 
for several models 0] . 

In the present paper we consider the shared informa- 
tion (or interdepency) in ground-states which arc super- 
positions of valence bond states. Our considerations ap- 
ply to ground-states in which the coefficients are all real 
nonnegative and therefore can be interpreted as probabil- 
ities of configurations [ij]. This implies that we consider 
the shared information resulting from correlations, in a 
bipartition of a classical and not a quantum system. The 
ground-states we study can describe equilibrium prob- 
lems (the spin 1/2, 577(2) symmetric one dimensional 



quantum chains [f| [f| are an example) but also proba- 
bility distribution functions (PDF) of stationary states 
of stochastic processes. We are going to concentrate on 
the latter and therefore also encounter systems which are 
scale invariant and not conformal invariant. We will show 
that if the system is conformal invariant, each entangle- 
ment estimator E(l,L) behaves like S q (l,L). The con- 
stant 7 has different values for different estimators. If 
the system is scale invariant but not conformal invariant, 
the Eq. © stays valid but the finite-size scaling function 
fl3J) is different: 

E(l,L)= 7 ln(Lg(l/L)), g(x)=g(l-x), (4) 

where g(x) ~ ax for small x and C = 7I11 a. Like the von 
Neumann entanglement entropy of quantum systems, the 
estimators detect the existence of long-range correlations. 

We present four estimators of the shared information 
and compare them considering two models defined using 
the same configuration space. We give here only the main 
results, all the details are going to be published elsewhere 
Q. Exact results are hard to obtain except for simple 
cases but for two estimators one can use Monte Carlo 
simulations for large system sizes and get reliable results. 
One of the estimators is not new fHH. {Vn. [Lsj . we are going 
to show its merits and limitations. 

In order to define the configuration space, we consider 
an open one-dimensional system with L sites (L even) 
connected by L/2 non-intersecting links (see Fig. |(TJ)). 
The links can be seen as U q {sl{2)) (a generalization of 
SU{2) [|) singlets. There are C L = L!/[(L/2)!(L/2+l)!] 
configurations of this kind. There is a bijection between 
link patterns and restricted solid-on-solid (RSOS) config- 
urations also called Dyck paths. A Dyck path is defined 
by taking L + 1 sites situated on the bonds of the link 
pattern. We attach to each site i non-negative integer 
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FIG. 1: (Color online) Example of a link pattern for L = 14 
and the corresponding Dyck path. In the latter there are four 
contact points and three clusters. The shared information is 
the largest in the left most cluster. 



heights hi which obey RSOS rules: 



H+i 



±1, h Q = h L = O' = 0, 1,...,L). (5) 



The height hi represents the number of crossed links at 
the site % (see Fig. [IJ. If hj = 0, at the site j one has 
a contact point. Between two consecutive contact points 
one has a cluster. There are four contact points and 
three clusters in Fig. [TJ It is easy to see that for a bipar- 
tition, large entanglements take place in large clusters. 
We present two models. In each of these models one has 
different probabilities for the various Dyck paths. In the 
cases in which one considers stationary states of stochas- 
tic models, the Dyck paths can be seen as an interface be- 
tween a substrate (h 2 i = 0, h-zi-i = 1, i = 0, 1, . . . , L/2)) 
covered by tiles (tilted squares) as shown in Fig. [T] and 
a gas of tiles (not shown in the figure) @. The PDF of 
the various Dyck paths are determined by the stochas- 
tic process. The latter is defined, in the time-continuous 
limit, by a non-hcrmitian Hamiltonian. 
We consider the following two models: 

A) The directed polymer model (DPM) (lol | . A configu- 
ration with m contact points gets a factor K m (K > 0). 
For K = 1, all configurations have the same probabil- 
ity 1/Cy, and represent the stationary PDF of the Rouse 
model [lH of a fluctuating interface. Using reflections 
about the horizontal axis, one can map the interface onto 
a random walker problem. The walker starts at the origin 
and crosses the horizontal axis after L steps. The den- 
sity of clusters is related to the first passage time problem 
and vanishes like L~ 3 / 2 for large L. In the whole domain 
< K < 2, one is in the same universality class as for 
K = 1. For K = 2 one gets a surface phase transition 
and for K > 2 the density of clusters stays finite in the 
thcrmodynamical limit. 

B) The stationary states of the raise and peel model 
(RPM). This is a stochastic model 0, [l2| in which the 
adsorption of tiles is local but the desorption is non-local. 
The model has a free parameter u. Typical configurations 
in the stationary state are shown in Figs. 11 and 15 of 
Ref. [l2j]. If < u < 1, the correlation length is finite 
and one has finite densities of clusters (see Fig. 12 in 
[l2l|). If u = 1, the average density of clusters vanishes in 
the thcrmodynamical limit and the system is conformal 
invariant (the dynamic critical exponent z — 1). This 



property makes the model special. The valence bonds 
represent U q (sl{2)) singlets for q = exp(i7r/3). The PDF 
in the stationary state has also remarkable combinato- 
rial properties. For u > 1 the system stays critical but 
conformal invariance is lost. The exponent z decreases 
smoothly with u from 1 to zero. There are fewer but 
larger clusters [l3[ than for u = 1. Because of its rich 
phase diagram, the RPM is an ideal playground to test 
various estimators. 

A bipartition of the system is obtained in the following 
way. The ensemble of Dyck paths (system C of size L 
(i + l sites)) is divided into two parts: the sites < 
i < I (part .4) and the sites I < j < L (part B). This 
implies the splitting of each Dyck path which at the site 
I has the height hi into two ballot paths One RSOS 
path which starts at i = and ends at the site / at the 
height hi and another one which starts at i = I, with 
height hi, and ends at a height zero at i = L. We denote 
by P(a(hi),b(hi)) the probability to have a given Dyck 
path in C formed by the ballot paths a(hi)(b(hi)) in A, 
respectively in B. We consider the marginals 



(6) 



and P(b(hi)). The probability to have a height hi at the 
site I is 



F(M = £iW0) = £W0)- 



(7) 



We present the four estimators. They all measure in 
different ways the amount of information that can be 
obtained about the ballot paths in B if one observes the 
ballot paths in A. 

I) Mutual Information. 



I(l,L) = - P(a(hi),b(hi)) In 

h t ,a(hi),b(hi) 

This is a known estimator (l5j . 
II) Boundary Shannon Entropy 



P{a{hi),b(h)) 
P(a(hi))P(b(hi))- 

(8) 



S{1, L) = H(L) - H(l) - H(L - I), 



(9) 



where H{M) = — J^k Pk ha(-Pfc) is the Shannon entropy 
for a system of size M and Pk is the probability to have a 
Dick path k. Notice that if, like in model A with K — 1, 
all configurations have the same probabilities and their 
number is Z(M), 

S(l,L) = -]nQ(L,l), 

Q(L, I) = Z(l)Z(L - l)/Z(L) - F{hi = 0), (10) 

where Q(L, I) is the probability to have the two systems 
A and B separated by a contact point at the site I. 
Ill) Density of Contact Points Estimator 



D(l,L)=ln(l/F(hi=0)). 



Ill) 
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Notice that D(l,L) and S(l,L) coincide if all configura- 
tions have the same probabilities. The physical meaning 
of D(l, L) is simple: if A and B have a small probabil- 
ity to be separated by a contact point the estimator is 
large. This should be the case since the shared informa- 
tion among A and B is large. 

In the continuum, F(hi — 0) can be replaced by p(l, L), 
the local density of contact points at the distance I from 
the origin for a system of size L. This is an average of 
a local operator. Let us observe that for 1 << I, L the 
density p stays finite, and therefore D(l,L) is also finite. 
If for large values of / and L one has, 



p(l,L) = l/{Lg(l/L)\ 



x 



(12) 



with g(x) ~ ax for small x (a being a constant), one 
obtains © and gj) with 7 = X. [The average number of 
clusters is p(L) ~ L~ x .] 

IV) Valence Bond Entanglement Entropy. 

This estimator was introduced independently by Chha- 
jlany et al [lj| and Alet et al [ijj and further studied by 
Jacobsen and Saleur [l8| (see also [HI). The estimator 
is the average height at the site I, for a system of size L: 



h{l,L)=Y,hiF{h). 



(13) 



We give the main results for the four estimators for 
each of the two models presented above (see Q). 
Directed Polymer Model: 
It is easy to show that for K = 1 one has 

1(1, L) = l/21n(Z(L-Z)/£)+7e-l/2(m(7r/2) + l) = 
~ 1/2 ln(/) + 0.303007 (L»l), (14) 

0. 



S(l, L) = D(l,L 



3 ML — 
2 ( T 



1 7T 

2 ln 8' 



(15) 



In (Tj"4")) 7 e is the Euler constant. Notice that Eq.© stays 
valid, the finite-size scaling functions in (T4")) and (jTSJ) 
are the same but different from the one given by ([3]) . We 
have checked Q that except for the additive constants, 
Eqs. (fT4"| and (TT5"]) are valid in the whole interval < 
A < 2 as expected from universality. The valence bond 
entanglement entropy does not get logarithmic but power 
corrections to the area law since h[l,L) ~ / 1 / 2 /(Z/L) in 
the whole interval < K < 2 [l0|. For A' > 2 and 
large values of L, the density of clusters is finite and all 
estimators verify the area law. 

Raise and Peel Model: 
u < 1 : 

D(l,L) and h(l,L) are finite since the average density 
of clusters is finite. 1(1, L) and S(l,L) were not com- 
puted. 

u=l (conformal invariance): 

1(1, L) and S(l,L) are hard to obtain since the PDF 
is know exactly only for small lattices. Some rough es- 
timates given in [?J show that they are compatible with 
©• 
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FIG. 2: (Color online). Left plot: The estimator h(l,L) for 
u = 1 for various lattice sizes L (I « L) as a function of 
InL Right plot: The scaling function [h(L/2, L) — h(l,L)] 
divided by ln(sin(7rZ/L)) for various values of u measured on 
a lattice of size L = 2400. For it = 1 one should get a constant 
independent on I. 



D(l , L) is obtained in the following way. As shown 
in (201 in an "almost" rigorous way, the density of con- 
tact points is the average of a local operator of a con- 
formal field theory and has the expression: p(l,L) = 
m/[L sin(7r//L)] 1 / 3 where m = -a/3F(-1/6)/(6tt 5 / 6 ). 
Using |JI2J| one finds: D(l,L) = l/31n(Lsin(7rZ/L)/7r) + 
0.28349. This is precisely Eq. © in which 7 = 1/3 is 
not given by the central charge of the Virasoro algebra 
(one would expect 7 = 1/6 if this would have been the 
case) but by the scaling dimensions of a local operator. 
Moreover one can estimate what would happen if the 
segment A would be inside an infinite system (two sep- 
aration points). The estimator D(l, 00) is given by the 
two-point correlation function of the densities separated 
by I, measured in fl~3j ] . One obtains a violation of the 
area law @ with 7 = 2/3. 

h(l, L) was obtained using Monte Carlo simulations for 
lattices up to L = 9600. The results are compatible with 
Eq. ©: h(l,L) = 0.277 ln(L sin(7r//L)/7r) + 0.73. These 
results were obtained in the following way. First we have 
taken I « L and plotted h(l,L) as a function of In I 
for various values of L (see Fig. [2]). One can see that 
there is a domain where we have a straight line which 
is L independent. This has allowed to get 7 and C (see 
Eq. ©). 

Next, we have considered the quantity [h(L/2,L) — 
h(l, A)]/ln[sin(7rZ/A)] for various values of L. If Eq. 
is valid one should obtain a constant equal to —7. The 
data are presented in Fig.[5]for L — 2400 and one can see 
that this is indeed the case. We should mention that con- 
sidering periodic boundary conditions and the boundary 
Coulomb gas formalism, Jacobsen and Saleur [l8| have 
obtained the value of 7 (one has to take half of their value 
since we deal with an open system) 7 = ^/S/2tt ~ 0.275 
which is compatible with our result. It is remarkable that 
the two estimators D(l, L) and h(l, L) give values for 7 
which are close to each other, 
u > 1 : 

The estimator D(l,L) can be computed using Monte 
Carlo simulations. A rough estimate of 7 can be obtained 
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using the equality 7 = X where the exponent X is re- 
lated to the density of clusters (see @). The exponent 
X varies between 1/3 and 1 when u increases from 1 to 
large values. One has X = 0.5, 0.6 and 0.85 for u = 1.2, 
1.5 and 10, respectively [l3l |. 

We have done a more detailed study for u = 4 (z ~ 
0.3) in this case and for D(l, L) we found (Ti"4"l) with 7 = 
0.73 ± 0.03 and a scaling function g{l/L) different of ([3]) 
(conformal invariance is lost at u = 4). We have also 
studied h(l,L) and found 7 = 0.63 ± 0.03 and a function 
g(x, L) (see right plot of Fig. F2]) equal within errors to 
the one observed for D(l,L) [21(. Notice that for both 
estimators the values of 7 have increased by more than 
a factor of two as compared with the values observed 
at u = 1. An increase of the shared information was 
expected since there are larger clusters connecting the 
subsystems A and B. 

The estimators defined above can be used for station- 
ary states of other processes not taking place in the Dyck 
paths configuration space. A simple example is the asym- 
metric exclusion problem (ASEP) with a density r of par- 
ticles on a ring of perimeter L. The role of the heights in 
the Dyck paths is played by the deviation of the number 
of particles in a subsystem (size I) from the number rl. 
The system being critical, one expects corrections to the 
area law. One finds indeed: 

1(1, L) = S(l,L) = D(l,L) = 
l/21n(Z(Z - l)/L) + l/21n(27rr(l - r)). 



Notice that 7 in (J2J) is r independent. 

We have shown that the four estimators of the shared 
information between two subsystems, defined above, ver- 
ify the area law. If the system is gapless, one obtains 
(with one exception) logarithmic corrections with the co- 
efficient 7 in ([2|) increasing if the shared information is 
larger. The exception is the average height, which in the 
directed polymer model gets power corrections probably 
due to the fact that the density of clusters decreases very 
fast with the size of the system. As a result, the existence 
of corrections to the area law can be used to detect the 
existence of phase transitions. Moreover, the observation 
of a finite-size scaling law like ([3|) can be an indication 
of conformal invariance. The estimators presented here 
have been generalized to the multi-partition case (see [7]). 
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